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We perform a detailed study of the modified gravity f{R) models in the light of the basic geo- 
metrical symmetries, namely Lie and Noether point symmetries, which serve to illustrate the phe- 
nomenological viability of the modified gravity paradigm as a serious alternative to the traditional 
scalar field approaches. In particular, we utilize a model-independent selection rule based on first 
integrals, due to Noether symmetries of the equations of motion, in order to identify the viability 
of f{R) models in the context of flat FLRW cosmologies. The Lie/Noether point symmetries are 
computed for six modified gravity models that include also a cold dark matter component. As it is 
expected, we confirm that all the proposed modified gravity models admit the trivial first integral 
namely energy conservation. We find that only the f{R) ~ (R!' ~ 2KY model, which generalizes the 
concordance A cosmology, accommodates extra Lie/Noether point symmetries. For this f{R) model 
the existence of non-trivial Noether (first) integrals can be used to determine the integrability of the 
model. Indeed within this context we solve the problem analytically and thus we provide for the 
first time the evolution of the main cosmological functions such as the scale factor of the universe 
and the Hubble expansion rate. 
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1. INTRODUCTION 

The comprehensive study carried out in recent years 
by the cosmologists has converged towards a cosmic ex- 
pansion history that involves a spatially flat geometry 
and a recent accelerating expansion of the universe (see 
[ll-Q and references therein). From a theoretical point of 
view, an easy way to explain this expansion is to consider 
an additional energy component with negative pressure, 
usually called dark energy, that dominates the universe at 
late times. In spite of that, the absence of a fundamental 
physical theory, regarding the mechanism inducing the 
cosmic acceleration, has given rise to a plethora of alter- 
native cosmological scenarios. Most of them are based ei- 
ther on the existence of new fields in nature (dark energy) 
or in some modification of Einstein's general relativity 
(GR) , with the present accelerating stage appearing as a 
sort of geometric effect ("geometrical" dark energy). 

The necessity to preserve Einstein's equations, inspired 
cosmologists to conservatively invoke the simplest avail- 
able hypothesis, namely, a cosmological constant, A (see 
for reviews). Indeed the so called spatially flat con- 
cordance ACDM model, which includes cold dark matter 
and a cosmological constant (A), fits accurately the cur- 
rent observational data and thus it is an excellent can- 
didate model of the observed universe. Nevertheless, the 
identification of A with the quantum vacuum has brought 
another problem which is: the estimate that the vacuum 
energy density should be 120 orders of magnitude larger 
than the measured A value. This is the "old" cosmo- 
logical constant problem (9|. The "new" problem 12] is 
related with the following question: why is the vacuum 



density so similar to the matter density at the present 
time? 

Such problems have inspired many authors to propose 
alternative dark energy candidates (see (l3| for review) 
such as A{t) cosmologies, quintessence. A:— essence, vector 
fields, pha ntom dark energy, tachyons and Chaplygin gas 
(see |9l. [l4| - [30y | and references therein). Naturally, in order 
to establish the evolution of the dark energy equation of 
state, a realistic form of H{a) is required which should be 
constrained through a combination of independent dark 
energy probes. 

On the other hand, there are other possibilities to ex- 
plain the present accelerating stage. For instance, one 
may consider that the dynamical effects attributed to 
dark energy can be resembled by the effects of a nonstan- 
dard gravity theory. In other words, the present acceler- 
ating stage of the universe can be driven only by cold dark 
matter, under a modification of the nature of gravity. 
Such a reduction of the so-called dark sector is naturally 
obtained in the f{R) gravity theories [3l|. In the original 
nonstandard gravity models, one modifies the Einstein- 
Hilbert action with a general function f{R) of the Ricci 
scalar R. The f{R) approach is a relative simple but still 
a fundamental tool used to explain the accelerated expan- 
sion of the universe. A pioneering fundamental approach 
was proposed long ago, where f{R) = R + mR"^ [3^ . 
Later on, the f{R) models were further explored from 
different points of view in [H, [s^ [s^ and indeed a large 
number of functional forms of f{R) gravity is currently 
available in the literature. It is interesting to mention 
here that subsequent investigations [35| confirmed that 
l/R gravity is an unacceptable model because it fails to 
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reproduce the correct cosmic expansion in the matter era. 

In this paper, we wish to test some basic functional 
forms of f{R) in the Ught of the Lie/Noether point sym- 
metries. The idea to use Noether symmetries in cosmo- 
logical studies is not new and indeed a lot of attention 
has been paid in the literature (see l37H47t ). Recently, 
we have proposed (see Basilakos et al. [4^) that the ex- 
istence of Lie/Noether point symmetries can be used as a 
selection criterion in order to distinguish the functional 
form of the potential energy V{(j)) of the dark energy 
models that adhere to general relativity (GR). In this 
work we would like to extend the paper of Basilakos et 
al. mi) by applying the same approach to f{R) models. 
In particular, the scope of the current article is (a) to 
investigate which of the available f{R) models admit ex- 
tra Lie and Noether point symmetries, and (b) for these 
models to solve the system of the resulting field equations 
and derive analytically (for the first time to our knowl- 
edge) the main cosmological functions (the scale factor, 
the Hubble expansion rate etc.). We would like to re- 
mind the reader that a fundamental approach to derive 
the Lie and Noether point symmetries for a given dy- 
namical problem living in a Riemannian space has been 
published recently by Tsamparlis & Paliathanasis [s^l (a 
similar analysis can be found in (49l - l55j ). 

The structure of the paper is as follows. The basic the- 
oretical elements of the problem are presented in section 
2, where we also introduce the basic FLRW cosmological 
equations in the framework of f{R) models. The geomet- 
rical Lie/Noether point symmetries and their connections 
to the f{R) models are discussed in section 3. In section 
4 we provide analytical solutions for those f{R) models 
which admit non trivial Lie/Noether point symmetries. 
Finally, we draw our main conclusions in section 5. 



2. COSMOLOGY WITH A MODIFIED 
GRAVITY 

Consider the modified Einstein-Hilbert action: 

1 



S- 



2fc2 



/ (R) + Cr^ 



(1) 



where £,„ is the Lagrangian of dust-like {pm = 0) matter 
and fc^ = SttG. Now varying the action with respect to 
the metric^ we arrive at 



= fc2 Tf 



2Df -{f-Rf 



(2) 



where the prime denotes derivative with respect to i?, 
is the Einstein tensor and T^f is the energy-momentum 



tensor of matter. Based on the matter era we treat 
the expanding universe as a perfect fluid which includes 
only cold dark matter with 4— velocity Ufj. . Thus the en- 
ergy momentum tensor becomes = —pm + (Pm + 
Pm)U^Uv, where pm and Pm = are the energy density 
and pressure of the cosmic fluid respectively. The Bianchi 
identity y** r^,y — leads to the matter conservation law: 



'iHpm = 



(3) 



the solution of which is pm = PmQO-~^ ■ Note that the 
over-dot denotes derivative with respect to the cosmic 
time t, o(t) is the scale factor and H = a/ a is the Hubble 
parameter. 

Now, in the context of a flat FLRW metric with Carte- 
sian coordinates 

ds^ = -df + a^{t){dx^ +dy^ + dz^) (4) 

the Einstein's tensor components are given by: 



Gl = -iH^, = -S^;: (2H + 3H 



(5) 



Inserting eqs.® into the modified Einstein's field equa- 
tions for comoving observers, we derive the modified 
Friedmann's equations 



3/ = fcVr, 



fR~f 



3Hf R 



(6) 



2f'H + 3f'H^ = -2Hf" R - (^f'" R^ + f" Rj - - — — 

^7) 

Also, the contraction of the Ricci tensor provides the 
Ricci scalar 

R^g^'R^.u^^i-^i.^^mH^ + H) . (8) 

Of course, if we consider /(i?) = R then the field equa- 
tions ^ boil down to the nominal Einstein's equations 
a solution of which is the Einstein de Sitter model. On 
the other hand, the concordance A cosmology is fully re- 
covered for f{R) = R - 2 A. 

From the current analysis it becomes clear that unlike 
the standard Fricdmann equations in Einstein's GR the 
modified equations of motion (|6]) and ^ are complicated 
and thus it is difficult to solve analytically. However, the 
existence of non-trivial Noether (first) integrals can be 
used to simplify the system of differential equations ^ 
and ([7]) as well as to determine the integrability of the 
system (see section 4). 



2.1. The f{R) functional forms 



^ We use the metric i.e. the Hilbert variational approach. 



In order to solve the system of eqs.® and ([7]) we need 
to know apriori the functional form of f{R). Due to the 
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absence of a physically well-motivated functional form for 
the f{R) parameter, there are many theoretical specula- 
tions in the literature. Bellow we briefly present various 
f{R) models whose free parameters, namely (m, n, Rc) > 
(0,0,0), can be constrained from the current cosmologi- 
cal data. 

• The power law model [H, [H, [s^l : 

f{R) = R-m/R" . (9) 

• The Amendola et al. '35^ modified gravity model: 



f{R)=R-mR,{R/R,)P 
with < p < 1. 
• The Hu & Sawicki jssl ] model: 

{R/Rc?" 



(10) 



f{R) = i? - mRc 



'{R/Rc)^^ + 1 

• The Starobinsky model: 

/(i?) = R- mRc 1 - (1 + R^/R^^ ~" 

• The Tsujikawa ^ model: 

f{R) = i? - TOi?ctanh(i?/i?c) . 



(11) 



(12) 



(13) 



• The generalization of the ACDM model (hereafter 
AbcCDM model 



f{R) = {R" - 2Ar 



(14) 



where the product be is of order of unity 0(1) and 
c > 1. The latter inequality is due to the existence 
of the matter epoch. 

Detailed analysis of these potentials exist in the lit- 
erature, includin g th eir confrontation with the observa- 
tional data (see [1J| for extensive reviews). We would 
like to stress here that within the context of the met- 
ric formalism the above f{R) cosmological models must 
obey simultaneously the some strong conditions (for an 
overall discussion see [ll]). Briefly these are: (i) /' > 
for i? > Rq > 0, where Rq is the Ricci scalar at the 
present time. If the final attractor is a de Sitter point 
we need to have / > for i? > _Ri > 0, where i?i is 
the Ricci scalar at the de Sitter point, (ii) / > for 
R>Rq>0, (hi) f{R) ^ R-2A for i? > i?o and finally 



(iv) < ^(r) < 1 at r = 



- -M. - 



f 



f 



Notice, that the power law f{R) model fails with re- 
spect to condition (ii) . The rest of the models satisfy all 
the above conditions and thus they provide predictions 
which are similar to those of the usual dark energy mod- 
els, as far as the cosmic history (presence of the matter 
era, stability of cosmological perturbations, stability of 
the late de Sitter point etc.) is concerned. Finally, in an 
appendix we discuss more f{R) models which however 
do not satisfy the conditions (i)-(iv) 62|. 



3. MODIFIED GRAVITY VERSUS 
SYMMETRIES 

In Basilakos et al. [1^ article we have proposed to use 
the Noether symmetry approach as a model-independent 
criterion, in order to classify the dark energy models that 
adhere to general relativity. The aim of this work is along 
the same lines, attempting to investigate the non-trivial 
Noether symmetries (first integrals of motion) by gen- 
eralizing the methodology of Basilakos et al. [4g] for 
modified gravity models (see section 2.1). This can help 
us to understand better the theoretical basis of the f{R) 
models as well as the variants from GR. 

In the last decade, a large number of experiments have 
been proposed in order to constrain dark energy and 
study its evolution. Naturally, in order to establish the 
evolution of the dark energy ("geometrical" in the cur- 
rent work) equation of state parameter a realistic form 
of H{a) is required while the included free parameters 
must be constrained through a combination of indepen- 
dent DE probes (for example SNIa, BAOs, CMB etc). 
However, a weak point here is the fact that the majority 
of the f{R) models appeared in the literature are plagued 
with no clear physical basis and/or many free parameters. 
Due to the large number of free parameters many such 
models could fit the data. The proposed additional crite- 
rion of Lie/Noether symmetry requirement is a physically 
meaning- full geometric ansatz, which could be employed 
in order to select amongst the set of viable models those 
which satisfy this constraint. Practically for those f{R) 
models which manage to survive from the comparison 
with the available cosmological data, our goal is to define 
a method that can further distinguish the f{R) models 
on a more fundamental (eg. geometrical) level and at the 
same time provides first integrals which can be used to 
integrate the modified Friedmann's equations. 

According to the theory of general relativity, the space- 
time symmetries (Killing and homothetic vectors) via the 
Einstein's field equations, are also symmetries of the en- 
ergy momentum tensor. Due to the fact that the f{R) 
models provide a natural generalization of GR one would 
expect that the theories of modified gravity must inherit 
the symmetries of the space-time as the usual gravity 
(GR) does. 

Furthermore, besides the geometric symmetries we 
have to consider the dynamical symmetries, which are the 
symmetries of the field equations (Lie symmetries). If the 
field equations are derived from a Lagrangian then there 
is a special class of Lie symmetries, the Noether symme- 
tries^, which lead to conserved currents or, equivalently, 
to first integrals of the equations of motion. The Noether 
integrals are used to reduce the order of the field equa- 
tions or even to solve them. Therefore a sound require- 



^ Note that the Noether symmetries are a sub-algebra of the alge- 
bra defined by the Lie symmetries [ssll . 
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ment, which is possible to be made in Lagrangian theo- 
ries is that they admit extra Noether symmetries. This 
assumption is model independent, because it is imposed 
after the field equations have been derived, therefore it 
does not lead to conflict with the geometric symmetries 
while, at the same time, serves the original purpose of a 
selection rule. Of course, it is possible that a different 
method could be assumed and select another subset of 
viable models. However, symmetry has always played a 
dominant role in Physics and this gives an aesthetic and 
a physical priority to our proposal. 

In the Lagrangian context, we can easily prove that 
the main field equations ([6|) and ([7]), described in section 
2, can be produced by the following Lagrangian: 



L = 6a/ + Qa^ f aR - 



(15) 



in the space of the variables {a, i?}. Using eq. p^ we 
obtain the Hamiltonian of the current dynamical system 



E = 6a/ 



E = 6a'^ 



6a^f aR 



fR-f 
6 



fR-f 



RHf 



(16) 



(17) 



Combining the first equation of motion (|6]) with eq.()17p 
we find 



E 



(18) 



The latter equation together with p„j = pmoo^ ^ implies 
that 



PmO 



E 
2P 



^7nPcr,0 — 



2fc2 



E = mmiil (19) 



where fi™ = Pmo/pcr,o, Pcr,o = SHg/k^ is the critical 
density at the present time and Hq is the Hubble con- 
stant. 

We note that the current Lagrangian eq. p^ is time 
independent implying that the dynamical system is au- 
tonomous hence the Hamiltonian E is conserved {dtE = 
^ = 0). Therefore, all the f{R) functions described in 
section 2.1 admit the trivial Noether symmetry, namely 
energy conservation as they should. 



3.1. Extra Lie and Noether symmetries 



proposed in Basilakos et al. [4^ article for the dark en- 
ergy models that adhere to GR. For the modified gravity, 
namely f{R) that lives into a 2-dimensional Riemannian 
space {a, R} and which is embedded in the space-time, 
how many (if any) of the previously presented functional 
forms (see section 2.1) can provide non trivial Noether 
symmetries (or first integrals of motion)? As an exam- 
ple, if we find a modified gravity model (or a family of 
models) for which its f{R) admits non-trivial first in- 
tegrals of motion with respect to the other f{R) cos- 
mological models, then obviously this model contains an 
extra geometrical feature. Therefore, we can use this ge- 
ometrical characteristic in order to classify this particular 
fjR) cosm olog ical model into a special category (see also 

[il EM, Hill). 

In order to compute the Lie/Noether point symme- 
tries of equations of motion (|6]) and ([7]) , we consider the 
Lagrangian'^ (jlSp as the sum of a kinetic energy and a 
conservative force field. The kinetic term defines a two 
dimensional metric in the space of {a,R}. Following 
standard lines (see [1^ and references therein) the two 
dimensional metric takes the form 



ds^ = Ylaf'dc? ^ Ua^f'da dR 
while the "potential" is 

V{a,R) = -a^f'R-f) . 



(20) 



(21) 



The signature of the metric eg. (1^01) is -f-l and the Ricci 
scalar is computed to be i? = 0, therefore the space is 
the 2-d Euclidean space^. Using the kinematic metric 
(PO)) we can utilize the plethora of results of Differential 
Geometry on coUineations to produce the solution of the 
Lie/Noether point symmetry problem. 

We recall that the special projective algebra of the of 
the Euclidean 2d metric ^IU\\ consists of the following 
vectors: 

= ada - 'i^dR , = -da ~ \^dR , = -^dR 



H^-l3a + ll,0R,A^^f'9a^lin: 



Be 



= ^Or , = ada , A* = ^Or 



a 

7' 

3 



= -a^fda + -a 
2 2 /" 



/" 



3 2/' 



where K are Killing vectors (K^'*^ are gradient), H is a 
gradient nomothetic vector, A are Affine coUineations 



Here we briefly present only the main points of the 
method used to constraint the f{R) models. In particu- 
lar, let us assume a modifled gravity f{R) cosmological 
model which accommodates a late time " accelerated" ex- 
pansion and it satisfles the strong conditions (i)-(iv) of 
section 2.1. We pose here a similar question with that 



^ In the appendix B we discuss the Noether symmetries in non flat 
/(i?) models. 

For the traditional dark energy models the signature of the two 
dimensional metric is —1 which means that the 2-d space is 
Minkowski H . Also all two dimensional Riemannian spaces 
are Einstein spaces implying that if ij = the space is flat. 
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and P are special projective collineations. These are ten 
vectors whereas the projective algebra of the two dimen- 
sional flat space consists of eight vectors [63|. It can be 
shown that the vectors K^,H^ are a linear combination 
of the affine vectors A^, / = 3, 4. 

Now we are looking for Noether symmetries beyond the 
standard one, dt- Utilizing the potential eq. (PT|) and the 
theorems 1 and 2 of 0, i55i| we find that among the f{R) 
models explored here (see section 2.1), only the Af,cCDM 
model [61j] with (&, c) = (1, |) admits extra Lie/Noether 
point symmetries. In particular the Lie point symmetries 



are 



Cie^* + C2e-^*)X2 (22) 



where the quantity Xj^^ is also Noether symmetry with 
gauge function 



9l.. 



= 9\/rn ^cie 



C2e 



a^jR - 2A 



where ci^2 are constants and m — 2A/3. If we relax 
the condition of c > 1 [6l| then we discover a second 
AfccCDM model with (6, c) = (1, |) that accommodates 
two Lie point symmetries the Xj^-^ and the Noether point 
symmetry 



Xr 



2^t _ _£2_ 2\/?Si 



cie2v^* + C2e 



m 



dt 



(23) 



with gauge function 



'^Vm {cie 



2,/mt 



C2e 



*)a^(i?-2A) s 



We have to mention here that the f{R) = {R - 2Ky/^ 
model does not satisfy the condition (ii), namely f"{R) > 
0. 

To conclude the discussion we would like to stress 
that the novelty in this work is the fact that among the 
current modified gravity models (see section 2.1) only 
the AfccCDM model |61[, which generalizes the concor- 
dance ACDM model, admits extra Lie/Noether point 
symmetries. This implies that the Af,cCDM model can 
be clearly distinguished from the other modified gravity 
models. Interestingly enough, the existence of the extra 
Lie/Noether point symmetries puts even further theoret- 
ical constrains on the free parameters of the AbcCDM 
model, (6,c) — (1, |) and {b,c) = (1,|)- From now on, 
we focus on the latter f{R) models and in the next section 
we provide for a first time (to our knowledge) analytical 
solutions. 



4. ANALYTICAL SOLUTIONS 

Using the Noether symmetries and the associated 
Noether integrals we solve analytically the differential 
eqs.® and ©. 



4.1. Ab.CDM model with (6,c) = (1, |) 



Inserting f{R) ^{R~ 2Kf/^ into eq. (fT5l) we obtain 



L = 9aVR-2Ad'^ + — hR+— (R + 4A) VR-2A 

2VR - 2A 2 

. (24) 

Changing now the variables from (a, R) to (x, y) via the 
relations: 

the Lagrangian and the Hamiltonian become 
L = xy + Vo (y^ + mxy) (25) 



E = ±y-Vo (y^ + mxy) 



(26) 



where Vq = 1/9 and rh = 6A. 

The equations of motion, using the the Euler-Lagrange 
equations, in the new coordinate system are 



X — 3Vby — mVox = 



y - mVoy = 0. 



(27) 



(28) 



The Noether point symmetries (1221) in the coordinate sys- 
tem {x, y} become 



X'l, = (cie'^*+C2e-"*)a^ 



(29) 



where lu = VmVo = ■\/2A/3 and the corresponding 
Noether Integrals are 



h = e^'y - we^'y 
/2 = e-'^^y + uje-^'^^y. 



(30) 
(31) 



From these we construct the time independent first inte- 
gral 



2„,2 



$ = /1/2 = y - y 



(32) 



The constants of integration are further constrained by 
the condition that at the singularity {t = 0), the scale 
factor has to be exactly zero, that is, a;(0) = 0. 
We consider the cases $ = and $ 7^ 0. 
A. Case = 0. 

We have the following sub-cases.: 
A.l. Ji = /2 = 

The solution of the system of equations ((27)) - (l28l) is: 



X (t) = xie"* -I- a;2e-"* , y (t) = (33) 

and the Hamiltonian constrain gives E = Q where 
xi^2 are constants. The singularity condition gives the 
constrain xi — —X2- At late enough times the scale fac- 
tor evolves as a^{t) (x x (t) oc a^ie"^*. However, this par- 
ticular solution is ruled out because it violates y{t) 7^ 0. 
A.2./1 = 0{l2^ 0) 
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The solution of the system ((271)- (HH) is: 



4.2. AbcCDM model with (6,c) = (1, |) 



P 



(34) 



(35) 



and the Hamiltonian constrain gives E — — ^2/2^ where 
x-j^ 2 Q-re constants. The singularity condition gives the 
constrain 



X 1 ~1~ Xo ^ 



4a; m 



= 



(36) 



At late times the solution becomes a^{t) cx x{t) 

4a; -^m 

A.3. /2 = (/i ^ 0) 

The solution of the system ((271)- ([Ml) is: 



4CJ^TO 



(38) 



and the Hamiltonian constrain gives E — xiIilu where 
xi^2 are constants. The singularity condition gives the 
constrain 



Xi + X2 + 



= 



This particular solution is not viable because in the mat- 
ter era we have e^'^* ~ implying that y{t) ~ 0. 
B. Case $ 7^ 

In this case the /i_2 ^ 0. The general solution of the 
system (gTll-dlll) is: 



h 



y{t) = ^e"* e 



x{t) = xice'^* + X2Ge~'^*+ 



1 



4m w 



(39) 



(40) 



The Hamiltonian constrain gives E ^ lo (xigIi ~ 2^20-^2) 
where x\q^2G are constants and the singularity condition 
results in the constrain 



1 



X\G + 2^20 + , _ n (-^1 + h) + = 0. 

Interestingly, one can show that the general solution in- 
cludes a proper matter era in which H{a) oc a"^/^ (see 
appendix C). Also, at late enough times the solution be- 

comes a^(i) cx x{t) cx ^^|^ e^"*. 



Despite the fact that the current f{R) model is phys- 
ically unacceptable due to f"{R) < 0, below we present 
its analytical solution for mathematical interest. In this 
case the Lagrangian eq. (Uni) of the f{R) = {R- 2 A) 
model is written as 



L = 



21a .0 21 . ^ 1 3 (i? - 16A) 



-aR a-^ 



(i?-2A)« 32(^_2A)^ 8 (i?_2A)« 

(41) 

We introduce the new coordinates (u, v) by means of the 
transformations: 



21 



and 



/^, i? = 2A+-g, (x,2/)^(0,0) 



1 u 



X — —;=UV , y — — = — . 

V2 V2v 



(37) In the coordinates (u, v) the Lagrangian is 



12 



L=-u^---v^ + Vo-u' + 2Vo— 42 
2 2 8 w' 

where fh = — 14A , K) = and the Hamiltonian: 



1.2 1" -2 



,12 



^=2" -2^^ -^°8" "2^- 



(43) 



The Euler-Lagrange equations provide the following 
equations of motion: 



U .2 VqUI 



y2 4 



u + AVq^ = 



,13 



2 1 , V- 
V + -iiv + 24Vo— r = 0. 

The Noether symmetries (l23l) become 

/Cl 2At '^2 2At 



(44) 
(45) 



(46) 



where A ~ ^y/mVo — ^yfA. The corresponding 
Noether Integrals are 



A 



(47) 



/_ = -e-^^'£; + e-^^*itM 4^ Ae^^^'u^ (48) 
A 

Following fs^l we construct the following time indepen- 
dent integral using a combination of the first integrals 
dH, dUl) and 



= WqV 



12 



(49) 
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The first integral (f> is called the Ermakov-Lewis invari- 
ant^. Using the Ermakov-Lewis Invariant, the Hamilto- 
nian (|43|) and equation (|44|) are written: 



1 



VqUI 



11 

2u2 



= E 



= 0. 



(50) 
(51) 



4 

The solution of (ISTI) has been given by Pinney [gJ] and 
it is the following: 



u{t) 



„2At 



U2e 



-2At 



where ui_3 are constants such as 



^ 4A2 (t 



(52) 



(53) 



From the Hamiltonian constrain ((50)) and the Noether 
Integrals (gll) , (gH) we find 

E = -2A-it3 , /+ ^ 2Xu2 , /- 2Aui. 

Replacing (|52p in the Ermakov-Lewis Invariant (UHl and 
assuming 7^ we find: 

« (t) = 2^0Ae-^W (4^0 + e-i2A(t) 

where 

A (t) = arctan 



■ 2A 

Then the solution is 



(uie 



2\t 



U3 



(54) 

AX^uiyQ^. (55) 



-2At 



U2e 4- 2U3) ' (56) 

where from the singularity condition a; (0) — we have 
the constrain ui + U2 + 2w3 = , or 

2E-{I++I_)=0. (57) 

At late enough time we find A{t) ~ Aq, which implies 
a'^{t) (xx{t) (x e^*. 

In the case where (j) = equations ([50]) . (ISTI) describe 
the hyperbolic oscillator and the solution is 

u{t) = sinhAt , 2£' = A^. (58) 

From the Ermakov-Lewis Invariant we have 



v{t) 



A sinh Xt 



^Xvi sinhAt- 12v/|K)|e"^^* / 
where vi is a constant. In this case the solution is 

- 

A sinh^ Xt \ 



1 



V2 \ Xvi sinh Xt - 12J\\^\e~ 



2Xt 



At late times the scale factor varies (t) cx x (t) 



(59) 



(60) 



^ An alternative way to compute the Ermakov-Lewis invariant is 
with the use of dynamical Noether symmetries [49l | . The corre- 
sponding dynamical Noether symmetry is Xjj = v?vdv 



5. CONCLUSIONS 

In the literature the functional forms of f{R) of the 
modified /(i?) gravity models are mainly defined on 
a phenomenological basis. In this article we use the 
Noether symmetry approach to constrain these mod- 
els with the aim to utilize the existence of non-trivial 
Noether symmetries as a selection criterion that can dis- 
tinguish the f{R) models on a more fundamental (eg. 
geometrical) level. Furthermore the resulting Noether 
integrals can be used to provide analytical solutions. 

In Basilakos et al. [i^ we have utilized the 
Noether symmetry approach to study the dark energy 
(quintessence or phantom) models within the context of 
scalar field FLRW cosmology. Overall the combination of 
the work of Basilakos et al. [48| with the current article 
provide a complete investigation of the Noether symme- 
try approach in cosmological studies. From both works 
it becomes clear that the Noether symmetry approach 
could provide an efficient way to discriminate either the 
"geometrical" (modified gravity) dark energy models or 
the dark energy models that adhere to general relativity. 
This is possible via the geometrical symmetries of the 
FLRW space-time in which both GR gravity and modi- 
fied gravity (or scalar field) live. 

In the context of f{R) models, following the general 
methodology of [SS^ (see also the references therein), the 
Noether symmetries are computed for 6 modified gravity 
models that contain also a dark matter component. The 
main results of the current paper can be summarized in 
the following statements (see sections 3 and 4): 

• We verified that all the f{R) models studied here, 
admit the trivial first integral, namely energy con- 
servation, as they should. 

• Among the 6 modified gravity models only the 
f{R) = {R'' - 2kY AbcCDM model with (6,c) = 
(1, |) provides a cosmic history which is similar to 
those of the usual dark energy models [see condi- 
tions (i)-(iv) in section 2.1] and at the same time 
it admits extra integrals of motion. In general, we 
propose that the f{R) models that simultaneously 
obey the conditions (i)-(iv), fit the cosmological 
data and admit extra Noether symmetries (integral 
of motions) should be preferred along the hierarchy 
of modified gravity models. Of course, one has to 
test the f{R) = (i?- 2A)^/2 model against the cos- 
mological data (SNIa, BAOs and CMB shift param- 
eter). Such an analysis is in progress and will be 
published elsewhere. Therefore the AbcCDM modi- 
fied gravity model appears to be a promising candi- 
date for describing the physical properties of "ge- 
ometrical" dark energy. We argue that although 
the AbcCDM model [61[ was phenomenologically 
selected in order to extend the concordance A cos- 
mology, it appears from the current analysis that 
it has a strong geometrical basis. 
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• Section 4 provides for a first time (to our knowl- 
edge) analytical solutions in the light of Af,cCDM 
model that include also a non-relativistic matter 
(cold dark matter) component. 
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Appendix A: Additional / (R) models which admit 
extra Lie/Noether point symmetries 

From the mathematical point of view and for the 
completeness of the present study, we would like to 
give the form of all / (i?) functions which admit extra 
Lie/Noether point symmetries but do not pass the con- 
ditions (i)-(iv). 

• If / (i?) is arbitrary we have the Lie point symme- 
tries 

ATli 2 ^ hdt + hada 

and the sole Noether point symmetry dt with 
Noether integral (constant of motion) the Hamil- 
tonian E. 

• If / (i?) ~ i?5 the dynamical system admits the 
extra Lie point symmetries 



1 



- -a 2 



da 



2R 



Or , Xl^ = t[-da 



1 



2R, 



Xl, = tdt - 2Rda 
and the extra Noether point symmetries 



Y 2i? 1 2R 

JiN2 — —da ^Or, Aat — t — Oa ~ t—^OR 

a a'' a 



Xn, = 2tdt + -ada' mdR. 
with corresponding Noether Integrals 

I,^^^{aVR) ,Is^t^^{aVR)-aVR 



h = 2tE ~ Ga^aVR - 6-^R. 

Vr 



• If / (i?) ~ i?8 the dynamical system admits the 
extra Lie point symmetries 

Xl, = 2tdt - AROr , Xl, = t^dt + t(^^da- 4i?9i^) 

and the extra Noether point symmetries 

Xn, = 2tdt + ^ da-mdR , Xn, = fdt+t da ~ mdR 

with corresponding Noether Integrals 



21 d / o 1 
H^2tE^--(a-R-^ 



21 d 



/g = t^E - —t— [ a^R-i ) + —a^R-t. 
8 at 



21 



• If / (R) ~ i?" (with n 7^ |, |) the dynamical sys- 
tem admits the extra Lie point symmetry 



Xls 



1 



-tdt 



1 



-Rdf 



2{n~l) ' n-1 
and the extra Noether point symmetry 



Xn, = 2tdt + ( -a (2n -I) da- ^.RBr 



with Noether Integral 

h 2tE-8na^R"-^a (2 - n)-'ina^R"-^R {2n - 1) (n - 1) 

Finally with the above analysis we would like to give 
the reader the opportunity to appreciate the fact that the 
Lie/Noether point symmetries provided in the current 
appendix, can be seen as an extension of those found by 
Vakili m. 



Appendix B: Noether symmetries in spatially non 
flat f{R) models 

In this appendix we study further the Noether sym- 
metries in non flat f{R) cosmological models. Briefly, in 
the context of a FRLW spacetime the Lagrangian of the 
overall dynamical problem and the Ricci scalar are 

L ^ Gf'aa^ + Gf'Ra^a + {f'R - /) - 6Kaf' 



a + K 



i? = 6 



where K is the spatial curvature. Note that the two di- 
mensional metric is given by eq. (l20l) while the " potential" 
in the Lagrangian takes the form V{a,R) = —a^{f'R — 
/) -I- Kaf. Based on the above equations and using the 
theoretical formulation presented in section 3, we find 
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that the f{R) models which admit non trivial Noether 
symmetries are: f{R) ^ (R - 2A)3/2, f[R) = i?3/2 
f{R) = R^. Notice, that the f{R) = {R - 2Ay/^ does 
not accept an analytical solution. 

In particular, inserting f{R) — {R—2A)^^^ into the La- 
grangian and changing the variables from (a, R) to (a;, y) 
[see section 4.1] we find 



L 



xy 



+ Vo {y^ + m.xy) - Ky 



E = xy — Vq (y + mxy) + Ky 



(Bl) 



(B2) 



Obviously, in order to reveal the evolution of the Hubble 
parameter in the matter era, in which the evolution of the 
matter density dominates the global dynamics, we have 
to understand the evolution of the first and the third term 
in eq. (ICip . Using R — 2K + ^ (see the transformations 
in section 4.1) we have, after some simple algebra, that 



■J/ 



-3 



9(i?- 2A)i/2 



y 

(C2) 



where K = 3{6^^^K). Therefore, the equations of motion 
are 

X - 3Foy^ - fhVox + K = 
y - fhVoy = . 

The constant term K appearing into the first equation of 
motion is not expected to affect the Noether symmetries 
(or the Integrals of motion). Indeed we find that the 
corresponding Noether symmetries coincide with those of 
the spatially flat f{R) = {R~2Af/^ model [see eqs. (|g^ . 
dSni), (ED), dSH)]- However, in the case of K ^ (or K ^ 
0) the analytical solution for the x-variable is written as 
XK{t) = x{t) + where x{t) is the solution of the flat 
model K = (see section 4.1). Note that the solution 
of the y- variable remains unaltered (see section 4.1 or 
equation C4). As expected, in the spatially flat regime 
K ^ the current equations reduce those equations of 
section 4.1. 

Similarly, in the case of /(i?) = R^/^ and f{R) = R^ 
the Noether symmetries can be found in appendix A. Of 
course, we again confirm that all the proposed modified 
gravity models with K accommodate the trivial first 
integral dtE — (energy conservation). 



f R 



R 



d{y^/x)/dt 



f 2(i?-2A) 2iyyx) 



(C3) 



For the benefit of the reader we repeat here the general 
solution of the system: 



' 2a; 2uj 



x{t) = xice"^* + X2Ge"'^*+ 
1 



(C4) 



(C5) 



where /1.2 7^ are the Noether Integrals. 

Inserting the general solution into eas. (|C2[) . (IC3|) and 
using at the same time that e""^* ~ we find 



,2Pm0a 2k -3/ ^ \1/2 



2k'^m 



Appendix C: Testing the analytical solutions 



In this appendix we would like to test the validity of 
our analytical solutions in the case of f{R) = {R—2Kf'/'^ 
model. Bellow we investigate the behavior of the Hub- 
ble parameter in the matter dominated era. First of all 
inserting pm = PmoOL~^ (see our eq fT5|) into the modified 
Friedmann equation (see eqlH) we get: 

,2Pmoa-' , f'R-f Hf"R 



3/' 



6/' 



/' 



(CI) 



/ R d{y^lx)/dt 



f 



2{y^/x) 



< 1 



Obviously, inserting the above results into the modified 
Friedmann equation ea. (|Cll) one can easily show that in 
the matter dominated era the evolution of the Hubble 
parameter tends to its nominal form namely H{a) — > 
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